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Abstract 

In India, six landscapes and source populations that are important for long-term conservation of Bengal tigers 
{Panthera tigris tigris) have been identified. Except for a few studies, nothing is known regarding the genetic 
structure and extent of gene flow among most of the tiger populations across India as the majority of them are 
small, fragmented and isolated. Thus, individual-based relationships are required to understand the species ecology 
and biology for planning effective conservation and genetics-based individual identification has been widely used. 
But this needs screening and describing characteristics of microsatellite loci from DNA from good-quality sources so 
that the required number of loci can be selected and the genotyping error rate minimized. In the studies so far 
conducted on the Bengal tiger, a very small number of loci (n = 35) have been tested with high-quality source of 
DNA, and information on locus-specific characteristics is lacking. The use of such characteristics has been strongly 
recommended in the literature to minimize the error rate and by the International Society for Forensic Genetics 
(ISFG) for forensic purposes. Therefore, we describe for the first time locus-specific genetic and genotyping profile 
characteristics, crucial for population genetic studies, using high-quality source of DNA of the Bengal tiger. We 
screened 39 heterologous microsatellite loci (Sumatran tiger, domestic cat, Asiatic lion and snow leopard) in captive 
individuals (n = 8), of which 21 loci are being reported for the first time in the Bengal tiger, providing an additional 
choice for selection. The mean relatedness coefficient (R = —0.143) indicates that the selected tigers were unrelated. 
Thirty-four loci were polymorphic, with the number of alleles ranging from 2 to 7 per locus, and the remaining five 
loci were monomorphic. Based on the PIC values (> 0.500), and other characteristics, we suggest that 16 loci (3 to 7 
alleles) be used for genetic and forensic study purposes. The probabilities of matching genotypes of unrelated 
individuals (3.692 x 10~ 19 ) and siblings (4003 x 10~ 6 ) are within the values needed for undertaking studies in 
population genetics, relatedness, sociobiology and forensics. 
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Background 

The conservation of the tiger, among the large felids, has 
been a global issue because of the extinction of three sub- 
species (Luo et al. 2004) and the decline of 93% of the 
habitat of the tiger (Karanth et al. 2010). The world tiger 
population is reported to have declined to as low a value as 
3200 (http:/ Avwf.panda.org/what_we_do/endangered_species/ 
tigers/tiger_resources/?188542/2010-Tiger-Factsheet) due to 
poaching as well as human activities that have resulted in 
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habitat fragmentation and depletion of wild prey species 
across the range of the species (Sunquist et al. 1999; Linkie 
et al. 2006; Sanderson et al. 2006). Among the different 
threats to the tiger, poaching and changes in landscape 
patterns are the greatest (Dinerstein et al. 2007; Goodrich 
et al. 2008; Walston et al. 2010), and hot spots of poaching 
may be identified by using genetic profile data, as has been 
done in tracking elephant ivory (Wasser et al. 2004). There- 
fore, a better understanding of the species at the individual 
level is needed for effective conservation planning and to 
avoid any further extinction of the extant sub-species. 

Among the extant tiger subspecies, the largest popula- 
tion (1706) is that of the Bengal tiger (Jhala et al. 2011), 
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which is the national animal of India and an endangered 
species listed under Schedule I of the Wildlife (Protec- 
tion) Act, 1972 of India. For science-based management 
of the species in fragmented landscapes, an understand- 
ing of its ecology, biology and genetics is crucial. The 
need for periodic monitoring of species genetics, espe- 
cially for large carnivores occupying highly exploited 
and fragmented landscapes, has also been emphasized 
(Anderson et al. 2004). Besides, reliable methods can 
be used to understand the causes responsible for the 
changing population demography are essential for de- 
signing the Tiger Conservation Plan (TCP) (Gopal et al. 
2007). However, with tigers being territorial, elusive, 
cryptic and nocturnal animals (Karanth et al. 2003), direct 
observation and enumeration are not feasible for obtaining 
such information. 

With the recent advances in molecular ecology, non- 
invasive genetic sampling and microsatellite markers for 
multi-locus genotyping have been used in studying ecol- 
ogy, biology and conservation genetics (Anderson et al. 
2004; Mondol et al. 2009a, 2009b; Miotto et al. 2007, 
2011; Charruau et al. 2011; Castilho et al. 2012; Reddy 
et al. 2012; Sharma et al. 2013), behavioral genetics 
(Langergraber and Vigilant 2011; Lyke et al. 2013) and 
socio-biology (Langergraber et al. 2013). 

Though, microsatellites have widely been used in un- 
derstanding genetics but a major constraint in the use of 
these loci is the need to isolate and characterize them 
using cloning and sequencing techniques. One of the 
ways of circumventing this step is to screen the varia- 
tions in microsatellites developed for other related spe- 
cies in order to find useful loci (Moore et al. 1991; 
FitzSimmons et al. 1995; Shepherd et al. 2002; Mantellatto 
et al. 2010). Therefore, numerous attempts have been 
made to use heterologous primers to support the conser- 
vation genetics of felids, viz. the jaguar {Panthera onca) 
(Ruiz-Garcia et al. 2006), snow leopard {Panthera uncia) 
(Waits et al. 2007), clouded leopard (Neofelis nebulosa) 
(Wilting et al. 2007), Siberian tiger (P. t. altaica) (Alasaad 
et al. 2011), cheetah {Acinonyx jubatus) (Charruau et al. 
2011), jaguarandi {Puma yagouaroundi) (Holbrook et al. 
2013), Indian leopard {Panthera pardus fused). (Mondol 
et al. 2009a; Dutta et al. 2012, 2013) and Bengal tiger 
(Bhagavatula and Singh 2006; Mondol et al. 2009b; Reddy 
et al. 2012; Sharma et al. 2013). However it is also useful 
to have large data available through screening of microsat- 
ellite loci across species. This will provide an alternate op- 
tion in selecting loci for a particular genetic study and 
may also lead to complement data or report if there are 
any discrepancies. 

Most of the studies undertaken so far on the Bengal 
tiger (Bhagavatula and Singh 2006; Mondol et al. 2009b; 
Reddy et al. 2012; Sharma et al. 2013) fail to provide de- 
tailed information on locus-specific genetic characteristics 



(polymorphic information content [PIC] and probability 
of identity [Pid]) and genotyping profile characteristics 
(stutter, allele to peak height etc.). Besides, information of 
these studies have been from fecal DNA, except for a few 
loci, which have been studied using high-quality DNA 
(Bhagavatula and Singh 2006; Mondol et al. 2009b). Thus, 
selection of the best loci for use in population genetics 
and forensic studies and minimizing genotyping errors has 
hitherto been precluded. 

Therefore, there is a strong need to describe locus- 
specific genotyping profile characteristics using DNA 
from a high-quality source, which has been suggested in 
the literature to minimize genotyping errors related to 
allele calling (Matsumoto et al. 2004). This has also been 
indicated in the guidelines of the ISFG (Gill et al. 2006, 
2012). Thus, we describe for the first time the screening 
and genotyping profile characteristics of 39 microsatel- 
lite markers developed for the Sumatran tiger {Panthera 
tigris sumatrae), domestic cat {Felis catus), Asiatic lion 
{Panthera leo persica) and snow leopard using DNA 
from a high-quality source. Of these, 21 loci are being 
reported for the first time in the literature for the Bengal 
tiger. Based on our findings, we suggest a combination 
of highly polymorphic dinucleotide and tetranucleotide 
repeat loci along with their genotyping profile character- 
istics for use in population genetic, forensic and non- 
invasive genetic sampling studies involving the Bengal 
tiger that will minimize allele calling errors by using 
locus-specific profile characteristics. Thus, the present 
study will provide better options in the selection and use 
of loci in population genetic and forensic studies carried 
out on Bengal tigers. 

Results and discussion 

Bengal tiger DNA samples (n = 8) were amplified suc- 
cessfully for all 39 heterologous loci, and data analysis 
using MICROCHECKER 2.2.3 (Van Oosterhout et al. 
2004) and GIMLET (Valiere 2002) clearly indicated the 
absence of null alleles, allele dropout, false alleles and 
scoring errors, associated with peak stuttering in geno- 
typing data. The mean value of the relatedness coeffi- 
cient (R = -0.143) also indicate that the selected tigers 
were not closely related to each other, as could be ex- 
pected in captive individuals. 

Three tetranucleotide repeat loci (Fca453, Fca731 and 
Fca749) and two dinucleotide repeat loci (6HDZ007 and 
Ple55) were found to be monomorphic in the Bengal 
tiger and were excluded from further analyses. In poly- 
morphic loci (n = 34), the observed allele size ranged 
from 78 to 315 bp (Table 1), whereas the number of al- 
leles (Na) per locus ranged from 2 to 7 (average 3.323). 
The effective number of alleles (Ne) per locus ranged 
from 1.438 to 4.923 (average 2.418). The average ob- 
served (H 0 ) and expected heterozygosities (H E ) for 



Table 1 Observed size range, genetic diversity statistics and genotyping profile characteristics for 39 microsatellite loci tested on 8 captive Bengal tigers 
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0.592 
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tetra-nucleotide 
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0.590 


0.539 


0475 
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Mean (based 
only on overall 
polymorphic 
loci, n = 34) 






3.323 


2.418 


0.625 


0.548 


0.482 


-0.103 
















Mean (based 






4.062 


3.081 


0.687 


0.664 


0.604 


0.022 

















only on 

suggested panel 
of polymorphic 
loci, n = 16) 

+ Locus recommended for panel of 16 microsatellite loci. 

Chr. Asn., chromosomal assignment of locus in species of origin; Nl, no information; T, tetranucleotide repeat; D, dinucleotide repeat; bp, base pairs; Na, number of alleles; Ne, number of effective alleles; H 0 , observed 
heterozygosity; H E , PIC, polymorphic information content; expected heterozygosity; P [D (locus), probability of identity between unrelated individuals; P !D Sibs (locus), probability of identity between siblings; Height ratio 
1, first stutter peak/main allele peak; Height ratio 2, minus A peak/main allele peak; Height ratio 3, plus A peak/main allele peak; Height ratio 4, heterozygote allele peak/main allele peak; PIC, polymorphic information 
content; F, Si inbreeding coefficients 'Williamson et al. (2002); 2 Menotti-Raymond et al. (1 999, 2005); 3 Singh et al. (2002); 4 Janecka et al. (2008); "significance of Hardy-Weinberg test (*P < 0.05). 
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polymorphic loci were 0.625 and 0.548, respectively. Four 
loci (PUN82, PUN100, PUN124, Ple57) had an H E level 
greater than 0.70. The higher value of H D compared with 
H E may be due to outbreeding that has probably taken 
place in a zoo as the animals were mixed from one popu- 
lation to another in India. A recent reduction in popula- 
tion size may cause a deficit of rare alleles compared with 
the number expected in a population at equilibrium. 
Since, rare alleles contribute comparatively little to H E , 
there will be an excess of H 0 while compared with a 
population at equilibrium among equal number of alleles 
(Cornuet and Luikart 1996; Garza and Williamson 2001). 
This hypothesis can be used to support the relatively high 
H 0 in this zoo population, in which the tigers may have 
different geographic origins. Hardy- Weinberg equilibrium 
(HWE) analysis may be affected by the small sample size, 
even though, we observed an HWE at all loci except 
6HDZ170, F85 and Fca723. Therefore, the higher values 
of heterozygosity are not due to HW disequilibrium. The 
method of Weir and Cockerham (1984) was used to calcu- 
late the inbreeding coefficient (F IS ), and the heterozygosity 
excess was examined at 68% of the loci (of which 74% 
and 26% were from dinucleotide and tetranucleotide 
repeat loci). The mean F IS value of the polymorphic 
loci was -0.103. Fifteen pairs of loci (6HDZ089 and F41, 
6HDZ089 and Fca272, F41 and Fca272, F42 and Fca304, 
F41 and Fca733, Fca272 and Fca733, Fca506 and Fca733, 
6HDZ089 and Fca740, F41 and Fca740, F53 and Fca740, 
Fca272 and Fca740, Fca506 and Ple57, Fca733 and Ple57, 
6HZ317 and PUN132, F41 and PUN132) indicated a sig- 
nificant linkage disequilibrium (LD) (P < 0.05). However, 
most of the selected loci of the domestic cat are located on 
different chromosomes. Therefore, the loci were selected 
carefully for the panel to avoid their physical linkage. The 
polymorphic microsatellite loci (n = 34) showed a mean 
PIC value 0.482, with sixteen loci having PIC values be- 
tween 0.511 and 0.770, six loci having PIC values between 
0.427 and 0.483 (which are considerably informative for 
population genetic studies (Botstein et al. 1980)) and the 
others having PIC values less than 0.400 (Table 1). 

The observed number of alleles indicates that the loci 
developed from the domestic cat, Asiatic lion and snow 
leopard have a greater number of alleles than do those 
from the Sumatran tiger (Figure 1). Pairwise statistical 
analysis (Mann- Whitney U test) indicates significant dif- 
ferences between the Sumatran tiger and domestic cat 
(P < 0.001), Sumatran tiger and Asiatic lion (P < 0.0001), 
domestic cat and snow leopard (P < 0.0001), Sumatran 
tiger and snow leopard (P < 0.0001 and Asiatic lion and 
snow leopard (P < 0.0001) but not between domestic cat 
and Asiatic lion (P < 0.105). This shows that the discrim- 
inatory power of the loci developed from the domestic 
cat, Asiatic lion and snow leopard is greater in Bengal 
tiger DNA samples. The majority of recent studies 



undertaken on felids have also used microsatellite loci 
developed for the domestic cat (Alasaad et al. 2011; 
Charruau et al. 2011; Dutta et al. 2012; Reddy et al. 
2012; Holbrook et al. 2013, Lyke et al. 2013; Sharma 
et al. 2013). Therefore, domestic cat microsatellite loci 
may enable a comparison of data across species to 
minimize ascertainment biases (Garner et al. 2005). 

The published reports indicate that there is a higher 
error rate for dinucleotide repeat loci than for tetra- 
nucleotide repeat loci during allele calling and this is dif- 
ficult to address due to a lack of genotyping profile 
characteristics (Cullingham et al. 2010). Therefore, we 
analyzed polymorphic dinucleotide repeat loci (n = 23) 
and tetranucleotide repeat loci (n = 11) separately to de- 
termine the level of allelic diversity, which has a strong 
significant role in individual identification. The number 
of alleles per locus at polymorphic dinucleotide repeat 
loci (n = 23) ranged from 2 to 7 (average 3.347), the aver- 
age observed and expected heterozygosities for these loci 
were 0.641 and 0.552, respectively, and the mean PIC 
value was 0.485 (Table 1). The number of alleles per 
locus at polymorphic tetranucleotide repeat loci (n = 11) 
ranged from 2 to 4 (average 3.272), the average observed 
and expected heterozygosities for these loci were 0.590 
and 0.539, respectively, and the mean PIC value was 
0.475 (Table 1). Our study clearly indicates that the poly- 
morphic dinucleotide and tetranucleotide repeat loci 
show more or less the same genetic diversity and other 
characteristics. Besides, there has been a choice of using 
tetranucleotide over dinucleotide loci to minimize prob- 
lems of allele calling (Cullingham et al. 2010). Thus, the 
domestic cat loci provide a better choice, with an ad- 
equate number of dinucleotide and tetranucleotide re- 
peat loci, compared with the loci developed for the tiger 
and other felids so far (Figure 1; Table 1). 

Allele scoring was easy for all the loci analyzed, and 
Figure 2 shows the allele scoring of one of the loci. 
Matsumoto et al. (2004) emphasized a need for inter- 
pretation of the locus-specific peak patterns and charac- 
teristics and suggested a novel algorithm for automated 
genotyping of microsatellites. We provide information 
for calculating the peak ratio of the first stutter, minus 
A, plus A and heterozygote allele (Table 1), which will 
make interpretation and allele scoring by others easier 
and more accurate. Such information are lacking for 
most of the studies so far undertaken for Bengal tigers. 

Hence, we suggest a panel of 16 microsatellite loci in- 
cluding polymorphic dinucleotide and tetranucleotide re- 
peat loci (Table 1) for genotyping-based studies carried 
out to understand the genetic structure of the population 
and to gather information on the ecology, biology and so- 
cial organization of the Bengal tiger from skin, tissue, fecal 
and hair samples. The suggested panel of 16 loci has 3 to 
7 alleles per locus (average 4.062); the average observed 
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and expected heterozygosities for these loci were 0.687 
and 0.664, respectively; and the mean PIC value was 0.604 
(0.511-0.770). Only two pairs of loci (F41 and PUN132, 
Fca506 and Ple57) showed a significant LD (P < 0.05), 
while chromosome location of PUN132 and Ple57 is not 



known (Table 1). Therefore, it should be checked whether 
they are also linked in other Bengal tiger populations. The 
mean F IS value of the suggested panel was also close to 
zero (0.022), which indicates that the selected captive 
population of Bengal tigers (n = 8) is in HWE. 
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The probability of identity (Pid), or probability of hav- 
ing the same genotype at multiple microsatellite loci of 
two individuals if they are drawn at random from a 
population, can be valuable information in a study where 
individual identification is needed. It can be estimated 
for differing number of loci (Waits et al. 2001). A P ID 
value of <0.01 (1 in 100) is considered essential for gen- 
etic studies in which population size estimation is re- 
quired (Mills et al. 2000). However, a sufficiently low P ro 
value of 0.001-0.0001 has been recommended in wildlife 
forensic applications for law enforcement (Waits et al. 
2001; Eiken et al. 2009; Lorenzini et al. 2011). A P ID level 
of <0.0001 has been used to study the population genetics 
of the bear and wolf (Waits et al. 2001). Figure 3 indicates 
that a combination of 5 polymorphic microsatellite loci 
from recommended panel (n = 16) was necessary to reach 
a Pid level of <0.0001 to adequately discriminate between 
individual tigers but was not sufficient for identifica- 
tion of siblings (Pid > 0.02). However, a combination of 
12-16 selected polymorphic heterologous microsatellite 
loci (Table 1) was adequate to reach a P ro level of <0.0001 
for discriminating siblings. The probability of identity of 
unrelated individuals determined using 16 polymorphic 
heterologous microsatellite loci was Pid (cumulative) = 3.692 x 
10~ 19 and of siblings P ro Sibs (cumulative) = 4.003 x 10" 6 , 
and thus it even meets the requirements of forensic studies, 
as suggested by Waits et al. (2001). The reported numbers 
of individuals in tiger populations in different protected 
areas of India range from 4 to 718 (Jhala et al. 2008), and 
some of the populations may be considered to be highly in- 
bred due to isolation and small population sizes. We rec- 
ommend the use of the suggested panel of 16 loci (Table 1) 
as it will not lead to any misidentification between two indi- 
viduals, including siblings, in small or inbred Bengal tiger 
populations. At the same time, a larger number of loci may 



introduce more genotyping errors when a low-quality 
source of DNA (viz. scat) is used (Creel et al. 2003). But 
the multiple-tube approach (Navidi et al. 1992; Goossens 
et al. 1998) and two-step multiplex PCR method can be 
used to overcome this problem without compromising 
the number of loci (Arandjelovic et al. 2009; Chang et al. 
2012), which are crucial for use in studies related to the 
ecology and biology of a species. 

When using different loci in studies involving samples 
that have been obtained non-invasively, the researcher is 
keen to know the error rates and amplification success 
rate. We tested the applicability of the recommended 
panel with noninvasive samples (scat) and blood from 
the same individuals and estimated the frequency of oc- 
currence of genotyping error rates. The values of the 
mean genotyping error rates were low and considerable 
for non-invasive genetic studies (allele dropout, 0.004 ± 
0.002 SD; false allele, 0.004 ± 0.002 SD and scoring error, 
0.006 ± 0.003 SD). These relatively low error rates may 
be due to the use of locus-specific profile characteristics, 
which leads to correct decisions in allele calling. We also 
did not observe any change or discrepancy in the genetic 
data compared with the data generated from blood 
samples. 

The key issue when using non-invasive genetic sam- 
ples, which are normally from poor-quality sources of 
DNA (especially scats), is identification and selection of 
loci that should have a higher amplification success rate as 
errors related to genotyping may be addressed by using 
other approaches that have been suggested (Matsumoto 
et al. 2004; Cullingham et al. 2010). We further tested our 
suggested panel of 16 markers and validated it with 50 
scat samples collected from different Bengal tiger popula- 
tions in India (Mishra et al. 2012). The preliminary results 
indicate that the average amplification success rate is 66% 




Figure 3 Probability of identity of unrelated individuals (P ID ) and probability of identity of siblings (P ID sibs) in locus combination 
using selected panel (n = 16). 
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with field-collected scat samples tested with a selected 
panel of 16 loci (Mishra et al. 2012), compared with 
other studies on carnivores, in which the reported success 
with fecal DNA is between 53% and 75% (Bellemain and 
Taberlet 2004; Bellemain et al. 2005; Smith et al. 2006; 
Murphy et al. 2007; Hansen et al. 2008). 

Our results of heterologous microsatellite loci, which 
have already been used in other studies, and additional 
loci (n = 21) will provide a wider choice for future efforts 
to assess the genetic diversity, existing range and genetic 
assignment of different populations of free-ranging 
Bengal tigers and minimize errors in allele calling. 

Material and methods 

Sample collection 

The first step before applying the non-invasive genotyping 
method to population monitoring and other aspects of the 
ecology and biology of the Bengal tiger is to identify a 
suite of hypervariable microsatellite loci using known 
good-quality tiger samples. To accomplish this, we ob- 
tained blood samples of 8 captive Bengal tigers which were 
sent to Wildlife Institute of India, Dehradun, India from 
Mahendra Chaudhury Zoological Park, Chhatbir, Mohali, 
India for DNA profiling. The histories of individual tigers 
and their translocation are inadequately documented in 
the Indian National Studbook for Bengal Tigers, 2011. 
Therefore, the place or geographic origin of these individ- 
uals is unknown. The reason behind opting for these indi- 
viduals in the present study is that if any microsatellite 
locus shows polymorphism in a captive population, that 
locus is supposed to show more polymorphism with wild 
individuals, which are thought to be outbred. DNA was 
extracted from their blood samples using Bio Robot EZ1 
(Qiagen, Germany). 

Scat samples from the same captive individuals (n = 8) 
and 50 scat samples from wild tigers were collected. A 
QIAamp DNA Stool Mini Kit (Qiagen, Germany) was 
used, following the manufacturer's protocol, to extract 
DNA from the scat samples. 

Selection, screening and genotyping of DNA from blood 
samples using heterologous microsatellite loci 

We selected and screened 25 dinucleotide and 14 tetranu- 
cleotide microsatellite loci that have been developed for 



the Sumatran tiger (Panthera tigris sumatrae) (Williamson 
et al. 2002), Asiatic lion (Singh et al. 2002), domestic cat 
(Menotti-Raymond et al. 1999, and 2005) and snow leop- 
ard (Janecka et al. 2008) to examine their allelic size range 
and polymorphism level in the Bengal tiger (Table 2). 
Polymerase chain reactions (PCR) were carried out 
in an Applied Biosystems 9700 thermocycler (Applied 
Biosystems, Germany) in a 10 pi reaction mixture 
containing 1 x PCR ABI Taq gold buffer, 2.0 mM MgCl 2 , 
0.4 mM dNTP mix, approximately 50 ng genomic DNA, 
4 pmol forward and reverse primers and 1 U Taq Gold 
DNA Polymerase (Applied Biosystems). Amplification was 
attempted for all 39 loci for all samples using PCR 
amplification conditions that have been published in the 
literature (Williamson et al. 2002; Singh et al. 2002; 
Menotti-Raymond et al. 1999, 2005; Janecka et al. 2008). 
The amplified PCR products were checked on 2% agarose 
gel in a 1 x TAE buffer. 

Statistical analyses 

The PCR products were scored on an ABI 3130 fluores- 
cence detection system using the GeneMapper software 
package (Applied Biosystems). The quality of the micro- 
satellite data was evaluated statistically for errors in geno- 
typing arising from null alleles (non-amplified alleles). 
Stutter peaks were scored using Micro-Checker 2.2.3 (Van 
Oosterhout et al. 2004). The frequencies of occurrence of 
large-allele dropout (short-allele dominance) and false al- 
lele were computed using GIMLET (Valiere 2002). To as- 
certain and obtain reliable genotypes, DNA from all eight 
blood, eight scats of captive Bengal tigers and fifty field 
collected scat samples were re-genotyped three to four 
times, respectively, at all the microsatellite loci screened 
so far (n = 39). Genetic diversity statistics for number of al- 
leles (Na), number of effective alleles (Ne), observed het- 
erozygosity (H 0 ) and expected heterozygosity (H E ) were 
generated using GenAlEx 6 (Peakall and Smouse 2006) 
and GENEPOP'007 (Rousset 2008). Using the allele fre- 
quencies, the polymorphic information content (PIC) 
of the markers was calculated using Cervus (ver. 3.0) 
(Kalinowski et al. 2007). The expected probability of 
matching genotypes for unrelated individuals (Pjx>) and 
siblings (Pid Sibs) was calculated for each locus using 



Table 2 Number of microsatellite loci, including dinucleotide and tetranucleotide repeats, from different species 
screened with Bengal tiger DNA samples 



S. No. 


Loci 
(n) 


Dinucleotide 


Loci repeats 

Tetranucleotide 


References 




Sumatran tiger 


8 


8 




Williamson et al. (2002) 




Asiatic lion 


4 


4 




Singh et al. (2002) 




Domestic cat 


20 


6 


14 


Menotti-Raymond et al. (19S 


9, 2005) 


Snow leopard 


7 


7 




Janecka et al. (2008) 
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GIMLET (Valiere 2002). GENEPOP'007 (Rousset 2008) 
was used to test the deviation from HWE. The F IS was 
determined using the probability test approach (Guo 
and Thompson 1992), with 10,000 dememorizations, 
500 batches and 10,000 iterations per batch in GENE- 
POP'007 (Rousset 2008). The inbreeding coefficients 
and the linkage disequilibrium (LD) were also tested 
using GENEPOP'007 (Rousset 2008). Considering the 
lack of details regarding individual tigers in the Indian 
National Studbook for Bengal Tigers, 2011, we estimated 
the Queller and Goodnight relatedness coefficients 
(Queller and Goodnight 1989) using GenAlEx 6 (Peakall 
and Smouse 2006). To ensure that the selected individ- 
uals were not related to each other, the level of relation- 
ship among the individuals was established using the 
R-value as suggested by Blouin (2003) and was calcu- 
lated using GenAlEx 6 (Peakall and Smouse 2006). 
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